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Updated predictions are presented for high energy neutrino and antineutrino charged and neutral 
current cross-sections within the conventional DGLAP formalism of NLO QCD using modern PDF 
fits. PDF uncertainties from model assumptions and parametrization bias are considered in addition 
to the experimental uncertainties. Particular attention is paid to assumptions and biases which 
could signal the need for extension of the conventional formalism to include effects such as ln(l/a;) 
resummation or non-linear effects of high gluon density. 
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I. INTRODUCTION 

(N 

Predictions of neutrino cross-sections at high energies have sizeable uncertainties which derive from the uncertainties 
on the parton distribution functions (PDFs) of the nucleon. In the framework of the quark-parton model, high energy 
neutrino deep inelastic scattering (DIS) accesses large values of Q 2 , the invariant mass of the exchanged vector boson, 
^v^j ' and small values of Bjorken x. the fraction of the momentum of the incoming nucleon taken by the struck quark. 
Thus in evaluating uncertainties on high energy neutrino DIS cross-sections it is important to use the most up-to-date 
information from the experiments at HERA, which have accessed the lowest x and highest Q 2 scales to date. HI and 
ZEUS have now combined the data collected in the years 1994-2000 to give very accurate inclusive cross-sections in 
the range 6 x 10~ 7 < x < 0.65 and 0.045 < Q 2 < 30000 GeV 2 [J. We consider modern PDF sets which include these 
Oh| data in order to provide the best benchmark cross-section for experiments searching for high energy cosmic neutrinos 
such as ANITA [|, IceCube H, the Pierre Auger Observatory [j] and RICE @, as well as forthcoming experiments 
rH ; such as ARA @, ARIANNA Q, JEM-EUSO ji and LUNASKA Q. 

Conventional PDF fits use the next-to-leading-order (NLO) DGLAP formalism [lol - fl3l | of QCD to make predictions 
• for DIS cross-sections of leptons on hadrons. At low x it may be necessary to go beyond the DGLAP formalism in order 
to sum ln(l/x) diagrams, as in the BFKL formalism fl4Hlal (fo r recent work see Refs. p7l42ll |). or to even consider 
non- linear terms as in the colour glass condensate model |22l |23|. The present paper is concerned with estimating the 

j» uncertainties on high energy neutrino cross-sections in the conventional NLO DGLAP formalism. However, the input 

. assumptions of some of the PDF fits are arguably beyond this formalism, and such cases will be highlighted. 

Pioneering estimates of high energy neutrino cross-sections were obtained at leading-order [24| using a PDF set 
(CTEQ4m) which is now well out of date, yet these values continue to be used in analysing current data (e.g. from 
ANITA 0). Two of the present authors (CSS) evaluated [26| the cross-sections at NLO using the ZEUS PDFs which 
included more modern HERA data [25[ and using a systematic procedure for estimating PDF uncertainties. More 
recently the cross-sections have been calculated in Ref. [27j (CTW) using the MSTW2008NLO PDFs, which use HERA 
data from the same vintage as the ZEUS PDFs and also account for PDF uncertainties. However, neither of these more 
recent works use the newly combined HERA results which have increased the accuracy of low x data by up to a factor 
of 3 (see Fig. 3 of Ref. [l|). It is the purpose of the present paper to re-evaluate the high energy cross-sections using 
the most up-to-date PDF sets, with particular emphasis on those which do use these precise, combined HERA data. 
The calculation is made using PDFs which were evaluated in NLO DGLAP fits, and our calculation of the neutrino 
structure functions and cross-sections is also made consistently at NLO. The input PDFs we consider treat heavy 
quarks by using general-mass- variable-flavour number schemes [28], [2|| , and we have also used such a scheme in our 
calculation of the structure functions. However, the difference between the use of a general-mass or a zero-mass scheme 
in the latter part of this procedure is negligible since the neutrino cross-sections are dominated by scattering at high 
Q 2 . We consider carefully the source of uncertainty on the input PDFs; these derive not only from the experimental 
uncertainty on the input data but also from model assumptions, and the form of the PDF parameterization. It is 
important to quantify these uncertainties carefully in order to be able to probe non-perturbative QCD effects at low 
x and/or new physics beyond the Standard Model [3ll - l34l | through a determination of the DIS cross-section using 
cosmic neutrinos which have energies extending up at least to ~ 10 11 GeV [35j |. 
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II. FORMALISM 



The kinematics of lepton hadron scattering is described in terms of the variables Q 2 , Bjorken x, and y which 
measures the energy transfer between the lepton and hadron systems. The double differential charged current (CC) 
cross-section for neutrino and antineutrino production on isoscalar nucleon targets is given by [36| 



d 2 a(v(9)N) G\Myj 



-a r (v(9)N), (1) 



dx dQ 2 4tt(Q 2 + M^) 2 x 
where the reduced cross-sections <j r (v{v)N) are 

a t (uN) = [Y + F»{x, Q 2 ) - y 2 F£(x, Q 2 ) + Y_xF£(x, Q 2 )] , (2) 



a T {yN) = [Y + F»(x, Q 2 ) - y 2 F£(x, Q 2 ) - Y_xF»{x, Q 2 )] , (3) 

and F2, xFs and Fl are related directly to quark momentum distributions, with Y± = 1 ± (1 — y) 2 . 

The QCD predictions for these structure functions are obtained by solving the DGLAP evolution equations at NLO 
in the MS scheme with the renormalisation and factorization scales both chosen to be Q 2 . These equations yield the 
PDFs at all values of Q 2 provided these distributions have been input as functions of x at some input scale Qq. 

In QCD at leading order, the structure function F^ is identically zero, and the structure functions F2 and xF% for 
charged current neutrino interactions on isoscalar targets can be identified with quark distributions as follows: 

F% = x(u + d + 2s + 2b + u + d + 2c), xFg = x(u + d + 2s + 2b - u - d- 2c), (4) 

and for antineutrino interactions, 

F% = x(u + d + 2c + u + d+2s + 2b), xF£ = x{u + d + 2c - u - d - 2s - 2b), (5) 

At NLO these expressions must be convoluted with appropriate co-efficient functions in order to obtain the structure 
functions (and F^ is no longer zero) but these expressions still give us a good idea of the dominant contributions. 
The contribution of the b quark will be suppressed until scales ~ m 2 , since the CKM element V t b ~ 1- Although 
the dominant contributions to the CC cross-sections come from Q 2 ~ m 2 , this does not mean that the b 

contribution is always suppressed, because the relevant scale for t production is the virtual boson-nucleon centre- 
of-mass energy, W 2 ~ Q 2 /x, and the high energy cross-sections are dominated by contributions from very small 
x ~ M 2 v /2m^E L , « M^/s (see Fig. [2]). This point had been overlooked in earlier work, including our own [261 ] . 

The neutral current cross-sections on isoscalar targets are given by expressions similar to Eqs. (JTJ [21 G2), with the 
W propagator replaced by the Z propagator, while the leading order expressions for the structure functions given by 

pr = * 

xFl' v = x[(u - u + d - d)(v u a u + v d a d )], 
where v u , Vd, a u , a d are the neutral current vector and axial- vector couplings for u— and d— type quarks. 



-(u + u + d + d) + (a 2 d + vj)(s + b + s + b) + (al + v 2 u )(c + c) 



(6) 



III. PARTON DENSITY FUNCTIONS 

The PDF4LHC group has recently benchmarked modern parton density functions 37]. Since our concern is with 
high energy neutrino cross-sections, rather than with LHC physics, we focus on PDF sets which make use of the 
newly combined accurate HERA data [lj]. Of all the PDFs considered by the PDF4LHC only HERAPDF1.0 [l|) 
and NNPDF2.0 [38[ used these data. However there has been a subsequent update of the CTEQ6. 6 |39| PDFs 
to CT10 jH which does use these data, while HERAPDF1.0 has recently updated to HERAPDF1.5 [4lfusing an 
preliminary combination of HERA data from 2003-2007 as well as the published combined data. We will utilise the 
CT10 and HERAPDF1.5 PDFs for the present study; we also consider the MSTW2008 PDFs in order to compare 
with other recent calculations of high energy neutrino cross-sections [27J, although we caution that these have not 
included the most accurate HERA low x data relevant to the present study. 

PDFs are generally determined by assuming a parameterization in x which is valid at a starting value of Q 2 = Qq, 
where the value of Qq is chosen to be sufficiently large that perturbative QCD calculations can be applied. The 
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form of the parameterization is usually Ax B (l — x) c P(x), where P(x) is some smooth function of x. Such forms are 
assumed for the light quarks and the gluon, whereas heavy quarks are generated dyamically from boson-gluon fusion. 
The PDFs at all other scales Q 2 > Qq are calculated by using the DGLAP equations to evolve the parameterized 
forms in Q 2 . The evolved PDFs are then convoluted with NLO matrix elements to calculate scattering cross-sections 
for processes of interest. These cross-section predictions are then fitted to deep inelastic data over a broad range of 
the x, Q 2 plane (HERA data cover 5 decades in both x and Q 2 ) in order to determine the parameters of the starting 
PDF parameterization. 1 In this way a small number of parameters (10-25, depending on the PDF set) are fitted 
to a large number (1000-2000) of data points. In this approach one can predict PDFs for x values below those for 
which data exists because the fitted functional form may be applied for all x, although uncertainties will naturally 
increase outside the fitted region. As we will illustrate, the uncertainty depends mainly on the theoretical prejudice 
underlying the parameterization at low x. 

The HERA data form the back-bone of PDF fits and are the only data which extend to the low x region. HERAPDFs 
use exclusively these data, while MSTW and CT also use older, fixed target data, Drell-Yan data including W and Z 
production, and Tevatron jet data. Moreover HERAPDFs are unique in using only proton data, so they are free of 
any assumptions concerning heavy target corrections and deuterium binding corrections [42j |. 

PDFs are presented accounting for the correlated systematic errors of the data as well as the statistical and 
uncorrelated sources, however, each PDF group has its own approach to the estimation of confidence limits on the 
uncertainties. A general discussion of approaches to estimating PDF uncertainties is given in Refs. [25|, [43| and the 
approaches used in the PDF sets considered here are reviewed in the PDF4LHC document [37[ • PDF uncertainties can 
also arise from input assumptions made in the PDF fitting. These include the form of the input PDF parameterization 
at the starting scale for evolution, Qg, the value of Q 2 itself, the kinematic cuts made on the data entering the fit, 
the value of a a (Mz), the values of the heavy quark masses — or even the scheme used to account for heavy quark 
production within the DGLAP fits. The PDFs considered here all account for the heavy quark production using 
general mass variable flavour number schemes, although the specific schemes differ: MSTW2008 and HERAPDF1.5 
use the Thorne-Roberts scheme [HHH, whereas CT10/CTEQ6.6 use the ACOT scheme 0. 

The experimental uncertainties on the PDFs are presented as eigenvector error sets. These eigenvectors represent 
linear combinations of the PDF parameters which are uncorrelated with each other since they are obtained by 
diagonalisation of the error matrix of the fitted parameters. Thus uncertainties on the PDFs, and quantities derived 
from them, can be calculated simply from adding in quadrature the difference between the independent eigenvector 
sets and the central PDF set. 2 However the confidence limits represented by these PDF uncertainties are not always 
set using the conventional tolerances of A\ 2 = 1 for a 68% c.L, and A\ 2 = 2.7 for a 90% c.l. Instead both CT(EQ) 
and MSTW increase the tolerances to account for the marginal inconsistencies of some of the input data sets, and for 
possible parameterization bias. The exact value of the tolerance is different for each eigenvector (see the individual 
PDF publications for details). The average tolerance for 90% c.l. is Ax 2 ~ 5 for MSTW2008 and Ax 2 ~ 10 for 
CTEQ6.6 and CT10. Note that MSTW supply both 68% and 90% c.l. sets, whereas CT10/CTEQ6.6 provide only 
90% c.l. sets; since we wish to consider 68% c.l. uncertainties, we simply scale the CT10 error sets by (1.64) -1 . 

The HERAPDFs use Ax 2 = 1 to set the size of their 68% c.l. experimental errors. The combination of the 
HERA data results in a consistent data set with small, well understood systematic errors 0. The HERAPDF1.0 
used NC and CC e + and er cross-section data from the first phase of HERA running 1992-2000 (HERA-I), while the 
HERAPDF1.5 update uses the preliminary combination of both HERA-I and HERA-II (2003-2007) running. Both 
of these HERAPDF sets provide not only PDF eigenvector error sets for experimental errors but also PDF variation 
sets which account for variation of input assumptions: namely variation of the charm and beauty masses, variation of 
the strangeness fraction in the sea, variation of the Q 2 cuts applied to the data, variation of the parametrization at 
Qq and of the value of Q\ itself. Not all of these sources of uncertainty are specifically considered by the other PDFs. 

A potentially important uncertainty is the variation of the minimum Q 2 cut for data entering the fit. This is set at 
Q 2 > 3.5 GeV 2 for HERAPDF 1.5 and at Q 2 > 2.5 GeV 2 for CT10 and MSTW2008. Although this cut should be high 
enough to avoid the non-perturbative region, it may include low x data for which BFKL, or other beyond-DGLAP 
effects, are already important. Hence we vary this cut to investigate the possible bias thus introduced. 

More importantly, HERAPDFs consider a variant of the input parameterization of the gluon at the starting scale 
which allows the gluon to become negative at low x, Q 2 . At NLO the gluon PDF does not have to be positive, 
although one might consider that it going negative signals a breakdown of the DGLAP formalism. This form of 
the gluon is in fact standard in the MSTW2008 parameterization; by contrast, the CT(EQ) analyses do not allow 



1 Note that not all parameters are fitted: the normalisations of the gluon and the u and d valence quarks are fixed by the momentum and 
number sum rules respectively while some other parameters are fixed by model assumptions; the individual publications of the PDF 
groups should be consulted for details [3?1| . 

2 We have used Eq. (19) of Ref. [45| . which allows for asymmetry of the positive and negative excursions along the same eigenvector 
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such negative gluon variants. However measurable quantities such as the longitudinal structure function Fl, which is 
closely related to the gluon at small x, must be positive. Therefore we check whether this is the case when a negative 
gluon PDF is encountered. 

A further uncertainty concerns the choice of the value of a s (Mz) because the shape of the low x gluon and the 
value of a s (Mz) are correlated through the DGLAP equations. All the PDF sets considered have been determined 
with a fixed value of a s (M z ) — 0.120 for MSTW2008, 0.118 for CT/CTEQ and 0.1176 for HERAPDF1.5. However 
PDF sets with a range of fixed values of a s (M z ) are also supplied for all these PDFs. The uncertainty on the PDFs 
due to a s (Mz) can be obtained by varying the value of a s (Mz) by ±0.0012 at 68%c.l. [46| . 

The PDFs discussed above are calculated at NLO in the DGLAP formalism and when they are used to make 
predictions of the neutrino DIS cross-sections the calculations should naturally be performed at a consistent order 
as we have done. We have also checked the impact of the use of NNLO PDFs and evolution for HERAPDF1.5 and 
found it to be negligible. 



IV. TECHNICAL DETAILS 

The calculation of the CC and NC cross-sections in NLO has been performed using DISPred 47|. QCDNUM [48[ as 
included in LHAPDF is used for the evaluation of the structure functions F2, Fl and xF%. The numerical integration 
of the differential cross-section is performed using the VEGAS Monte Carlo integration routine from the GSL library (4|| , 
reducing the integration error to less than 1%. 

The most up-to-date PDFs in the LHAPDF [5(| format are given in form of grids with a limited kinematic range. 
For example, the grid extends down to 10~ 8 in x for the CT10 and HERAPDF1.5 sets and down to only 10 -6 for the 
MSTW2008 set. Below this range the PDFs 'freeze', i.e. they take the value at the lower grid boundary. However, 
for very high energies the total cross-section has significant contributions from lower x (see also Fig. [3]) such that the 
cross-section calculated using the LHAPDF sets will be too small. Therefore, we have used implementations for the 
PDFs provided by the CTEQ and MSTW [52[ groups which allow an extrapolation beyond the grid. Of course, 
such an extrapolation using e.g. polynomials is not necessarily physical but in the absence of a grid that extends to 
low enough values of x or the parametric form of the PDFs, this is the only viable alternative. We have checked the 
agreement between the cross-sections calculated with 'freezing' and with the extrapolations and the difference is more 
than 50% for MSTW2008 but only 3% for CT10 at E v = 5 X 10 11 GeV. 



V. RESULTS 



Figure [T] shows the predicted sea and gluon distributions from the HERAPDF1.5 fit and their fractional uncer- 
tainties, at various Q 2 values. This shows the importance of the low x contribution and illustrates that the PDF 
uncertainties are largest at low Q 2 and at low x. PDF uncertainties are also large at very high x but this kinematic 
region is not important for scattering at high neutrino energies. 

This is seen from Figure [2] which indicates the kinematic regions in the x, Q 2 plane which contribute to the neutrino 
cross-sections for two representative values of the neutrino energy: s = 10 s GeV 2 (=>■ E v = 5.3 x 10 7 GeV) in the 
left panel and s — 10 10 GeV 2 E v = 5.3 x 10 9 GeV) in the right panel. (We do not show the antineutrino 
cross-sections separately because these are very close to the neutrino cross-sections at high energy.) One can see that 
the dominant contributions come from 500 < Q 2 < 50000 GeV 2 and 10 -6 < x < 10~ 2 for the lower neutrino energy, 
and 10 -8 < x < 10 -4 , for the higher neutrino energy. Thus although PDF uncertainties are large at low x and low 
Q 2 , we see that the dominant contributions to the cross-section do not come from the kinematic regions of greatest 
uncertainty for the PDFs. 

Figure [3] compares the gluon PDF and its uncertainty at Q 2 = 10000 GeV 2 for the three PDFs which we consider. 
This value of Q 2 is in the middle of the range which contributes significantly to the neutrino cross-sections. We see 
that the central values of the gluon PDFs are all very similar, whereas the uncertainty estimates differ. The CT10 and 
HERAPDF1.5 uncertainties are actually very similar if we leave out member 52 from the CT10 error set. This error 
set was introduced into the CT10 analysis to allow for a larger uncertainty at low x 53]. Previous CTEQ analyses 
such as CTEQ6.6 [39( do not have such an extreme error set — see left panel of Fig. 01 The problem with such an 
ad hoc introduction of a steeply increasing gluon PDF is that at low x it leads to a very strong rise of the neutrino 
cross-section which seems unphysical (see later discussion). 

The larger error band of MSTW2008 is partly due to the fact that it does not include the most up to date HERA 
data, which have significantly reduced errors at low x. However the more striking difference between MSTW2008 
and both HERAPDF1.5 and CT10 is the downward divergence of its error band which is due to the gluon becoming 
negative at low x, Q 2 . This is best understood by reference to the right panel of Fig|4]which shows HERAPDF1.5 both 
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FIG. 1: The HERAPDF1.5 gluon PDFs and their fractional uncertainties — from the experimental errors only (EIG), and 
from both experimental errors as well as model and parameter variations (EIG & VAR) — at various values of Q 2 . 
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FIG. 2: Kinematic regions of the x, Q plane and their contribution to the total neutrino cross-section using HERAPDF1.5 for 
two different values of s. The labels show the relative contribution to the total cross-section contained within each contour. 
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FIG. 3: Left panel: Gluon structure function at Q 2 — 10 4 GeV 2 for the three PDF sets used. Right panel. The relative 
deviations and uncertainties (at 68% c.l.) with respect to the central value of HERAPDF1.5. The uncertainty bands are shown 
with member 9 for HERAPDF1.5 and member 52 for CT10. 



with and without member 9 (the variant which allows the gluon to become negative at low x, Q 2 ). However Fig. [S] 
shows that while the longitudinal structure function Fl is always positive for HERAPDF1.5, it becomes negative, 
hence unphysical, for some of the MSTW2008 PDF error sets. 

The total neutrino cross-sections are now obtained by integrating the predicted double differential cross-section 
d 2 a/dxdy with no cuts on either kinematic variable. Fig. [5] and [7] show the NC and CC neutrino cross-sections as a 
function of E v as evaluated from the HERAPDF1.5 and CT10 PDFs respectively. The PDF uncertainties of these 
predictions are shown relative to the central values underneath each plot. 

It is clear already that the variation of input assumptions considered by the HERAPDF1.5 — apart from the negative 
gluon variation — do not increase the error band very significantly. Thus the other parameterization variations and 
the input values of the charm and beauty mass and the fraction of strangeness in the sea are not important for the 
predictions of neutrino cross-sections. However it is interesting to consider the variation due to the choice of the 
minimum Q 2 for data entering the fit. Decreasing this cut to Q 2 > 2.5 GeV 2 has negligible effect but the effect of 
increasing the cut to Q 2 > 5 GeV 2 is just about visible in Fig. \§\ This larger Q 2 cut also cuts out data at x < 10~ 4 
— the kinematic region where there are already some hints of beyond-DGLAP behaviour such as BFKL or non-linear 
effects. This results in a steeper low x gluon and we can see this in the marginally enhanced neutrino cross-section. 

Another potentially important effect comes from the variation of a s (Mz) which is correlated to the gluon PDF such 
that lower values of a s (Mz) result in a steeper low x gluon. We evaluate this by considering a 68% c.l. variation of 
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FIG. 5: Left panel: The longitudinal structure function F L calculated from member 18 of MSTW2008nlo68cl for s = 10 8 GeV 2 
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Aa a (Mz) ± 0.0012 from its central value. The slight enhancement in the neutrino DIS cross-section is so small that 
it is not noticeable on the scale of Figs. [5] and [7] hence we have not attempted to show it. 

In Figs[8](top panels) we compare the NC and CC cross-sections, along with their total uncertainties (including that 
coming from the variation of a s (M z )), as predicted by HERAPDF1.5 and CT10. The MSTW2008 central prediction 
is also included for comparison. In Fig [5] (bottom panels) we emphasize the small differences in the central values of 
the PDFs and their relative uncertainties. In order to highlight the effect of the extreme members of HERAPDF1.5 
and CT10 in Figs [9l we show these plots without member 9 of the HERAPDF15 variations (which allows for the 
gluon to become negative at low x and Q 2 ) and without member 52 for CT10 (the cross-section for which rises oc E®' 7 
whereas for the central member it rises oc E®' 3 ). However any power-law rise in the cross-section will eventually 
violate the Froissart bound, which requires the rise to be no faster than lo g 2 s [5~j ] . This should result in a reduction 
of the cross-section at high energies, by a factor of ~ 2 at E v — 10 12 GeV |55| and perhaps even more (56|. 

There are small differences from the previous work of CSS [26[ and CTW [27[. In Fig.[TU]we compare our calculation 
using MSTW2008 PDFs to that of CTW. We find that we can reproduce their results well only if we use the MSTW2008 
NLO PDFs together with an leading- order treatment of the coefficient functions, rather than with a consistent NLO 
approach. In Fig[Tl]we compare our calculation using HERAPDF1.5 to that of CSS [2(|. The predictions for both 
the central values and the uncertainties of the neutrino NC cross-section are quite close. This is also true for the CC 
cross-section at low E v , however above 10 4 GeV the difference increases since in the present work we have included 
the contribution of the b quark which was missed out in our earlier work. Accordingly the bounds derived on the 
cosmogenic neutrino flux using the CSS cross-sections 0, Q are conservative, being ~ 20% too high. 

In Fig . [T2] we compare our results for HERAPDF1.5 to the cross-sections used in the neutrino event generator 
ANIS [57] which is based on CTEQ5D. Note that at energies below a TeV (which is the most important energy range 
for neutrino telescopes like IceCube [11]) there is a ~ 10% discrepancy. We also compare the CC cross-section for 
HERAPDF1.5 to its value in the GENIE low energy neutrino event generator [59| at around 100 GeV, finding the 
match to be consistent within errors. 
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FIG. 6: Neutrino DIS cross-section for CC and NC scattering as predicted by the HERAPDF1.5. Total relative PDF uncer- 
tainties including experimental errors as well as model and parameter variations (EIG & VAR) are shown beneath each plot, 
both with and without member 9. The marginal effect of increasing the Q 2 cut from 2.5 to 5 GeV 2 is also shown (member 8). 




FIG. 7: Neutrino DIS cross-section for CC and NC scattering as predicted by the CT10 PDF. Relative PDF uncertainties from 
experimental errors only (EIG), and also including parameter uncertainties (EIG & VAR) are shown beneath each plot, both 
with and without member 52. 
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FIG. 8: Comparison of the total cross-section (top panels) and uncertainties (bottom panels) for CC and NC scattering as 
predicted by the HERAPDF1.5, CT10 and MSTW2008 (central member only) PDF sets. The cross-sections and deviations 
for member 9 of HERAPDF1.5 and member 52 of CT10 are indicated by the dashed and dot-dashed lines, respectively. 
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FIG. 9: Same as Fig. [5] but excluding member 9 of the HERAPDF1.5 set and member 52 of the CT10 set. 
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FIG. 10: Our results for the neutrino DIS CC and NC cross-section using the MSTW2008 central member, relative to the 
results of CTW [27|. The calculation is done consistenly in NLO as well as, for illustration, using a LO code. 
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FIG. 11: The relative deviation of the cross-section calculated by CSS [2(| from our result for the HERAPDF1.5 central 
member. For the CC cross-section (left) we compare to HERAPDF1.5 with (upper pa nel) and without (lower panel) the 
6-quark contribution. For NC scattering (right) the 6-quark was included already by CSS [2g| and the agreement is excellent. 



VI. CONCLUSIONS 



We find that the predictions of high energy neutrino DIS cross-sections from the central values of HERAPDF1.5, 
CT10 and MSTW2008 PDFs are very similar. However the predictions for the uncertainties (deriving from the 
uncertainties on the input PDFs) differ quite strongly. In fact PDF uncertainties derive from the input assumptions 
as well as from the input experimental data. If we exclude error sets which either lead to too steep a rise in the cross- 
section, or allow the low x gluon to be negative at low Q 2 , then we find that the uncertainty estimates of HERAPDF1.5 
and CT10 — both of which use the most up-to-date, accurate HERA data — are remarkably consistent. 

Our results for the high energy neutrino and antineutrino CC and NC DIS cross-sections and their uncertainties 
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FIG. 12: The relative deviation of the ANIS ^ and GENIE [m| cross-sections from the HERAPDF1.5 central member. 
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FIG. 13: Neutrino and anti-neutrino cross-sections on isoscalar targets for CC and NC scattering according to HERAPDF1.5. 



using HERAPDF1.5 at NLO are shown in Fig.Q2J The general trend of the uncertainties can be understood by noting 
that as one moves to higher neutrino energy one also moves to lower x where the PDF uncertainties are increasing. 
The PDF uncertainties are smallest at 10~ 2 < x < 10 _1 , corresponding to s ~ 10 5 GeV 2 . Moving to smaller neutrino 
energies brings us into the high x region where PDF uncertainties increase again. This effect is greater for the 
HERAPDF1.5 because the HERA data have less statistics at high x than the fixed target data which are included 
in CT10; however these data have further uncertainties that are not fully accounted for in CT10, e.g. heavy target 
corrections, deuterium corrections and assumptions regarding higher twist effects. When the high x region becomes 
important the neutrino and antineutrino cross-sections are different because the valence contribution to xF% is now 
significant. This is seen in Fig. [T31 as is the onset of the linear dependence of the cross-sections for s < M^. Note 
that our predictions are made for Q 2 > 1 GeV 2 since perturbative QCD cannot sensibly be used at lower values. 
Moreover for s below ~ 100 GeV 2 , there can be contributions to the cross-section of 0(10%) from even lower values 
of Q 2 which are not accounted for here; hence we do not show results for E v below 50 GeV where there are other 
contributions to the neutrino cross-section and the use of a code such as GENIE 59] is appropriate. For higher 
energies, we intend to upgrade ANIS [57j to use the HERAPDF1.5 (differential) cross-sections. Meanwhile we have 
provided the total DIS cross-sections for CC and NC scattering of neutrinos and antineutrinos on isoscalar targets 
in Tables U and [TT] and recommend these as a benchmark for use by experimentalists. These cross-sections as well 
as those for isoscalar targets are available from a webpage [fijj; differential cross sections are available upon request. 
Any measured deviation from these values would signal the need for new physics beyond the DGLAP formalism. 
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TABLE I: Neutrino CC and NC cross-sections on isoscalar targets, along with their uncertainties, in the perturbative DGLAP 
formalism at NLO, using HERAPDF1.5 (both with and without member 9). 
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TABLE II: Antineutrino CC and NC cross-sections on isoscalar targets, along with their uncertainties, in the perturbative 
DGLAP formalism at NLO, using HERAPDF1.5 (both with and without member 9). 
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